3D MHD Simulations of accreting neutron stars: 
evidence of QPO emission from the surface. 
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Abstract. 3D Magnetohydrodynamic simulations show that when matter accretes onto neutron 
stars, in particular if the misalignment angle is small, it does not constantly fall at a fixed spot. 
Instead, the location at which matter reaches the star moves. These moving hot spots can be 
produced both during stable accretion, where matter falls near the magnetic poles of the star, and 
unstable accretion, characterized by the presence of several tongues of matter which fall on the 
star near the equator, due to Rayleigh-Taylor instabilities. Precise modeling with Monte Carlo 
simulations shows that those movements could be observed as high frequency Quasi Periodic 
Oscillations. We performed a number of new simulation runs with a much wider set of parameters, 
focusing on neutron stars with a small misalignment angle. In most cases we observe oscillations 
whose frequency is correlated with the mass accretion rate M. Moreover, in some cases double 
QPOs appear, each of them showing the same correlation with M. 
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MOVING HOT SPOTS AND QPOS 

In the last few years 3D MHD simulations have shown that in many cases, in particular 
if the misalignment angle is small (< 15°), matter falls in a ring around the magnetic 
pole. The brightest spot (the hot spot) of the ring is not fixed, but moves with an angular 
velocity related to the Keplerian velocity in a zone of the disk outside the magnetospheric 
radius. This is observed both in the case of stable accretion (e.g. J3) and in the case of 
unstable accretion (e.g. 0,fl3)). During unstable accretion the hot spots produced by 
instabilities in the equatorial region move, with a velocity similar to that of the inner 
disk (E^GO). The motion in this case is less regular and shorter. We studied the effects 
of these movements in detail with Monte Carlo techniques and showed that the features 
produced by those movements are quasi periodic oscillations. 

Correlations with the mass accretion rate 

We performed a number of simulations in which the only parameter changed was the 
mass accretion rate. All the systems reproduced in the simulations have M — 1.4M , 
R = 10km and misaligment angle = 2°. One system has a rotational period T = 4. 1ms, 
the other t = 5.4ms. In Fig. [T] we see the rotational frequency of the hot spots plotted 
against the mass accretion rate M. In one case the accretion is only on the poles, in the 
other both in the poles and in the equatorial zones due to instabilities. In the latter case, 
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FIGURE 1. Frequency variation of the hot spot movements with the mass accretion rate |4|. Results 
are scaled for the two cases B = 10 8 G and B = 5 • 10 8 G. The first case (right) shows stable accretion from 
the magnetic pole, with matter falling onto the surface along the funnel formed by the magnetic field 
lines. The hot spot moves around the magnetic pole with a certain angular velocity, which we plot against 
the mass accretion rate M, finding a good correlation. The second case (left) presents instabilities at the 
magnetospheric boundary (or truncation radius). In this case we have both accretion onto the magnetic 
pole with the usual funnel-like stream, that moves as in the first case, and accretion from tongues that 
pierce the magnetosphere falling near the equator, moving approximately at the frequency of the inner 
disk. Hence the "double track" in the plot. 



the hot spots appearing in the two accretion zones move together with the accretion rate. 
This result is, of course, extremely interesting to interpret the presence of multiple QPOs 
in accreting neutron stars. We have a mechanism to produce very bright pairs of QPOs, 
whose spectrum can easily be different from the background emission from the disk, and 
whose characteristic frequencies change together with M. Moreover, the quality factors 
of the two peaks show a different behavior, as observed in many works (e.g. (6]],0). 
Submitted to MNRAS. A preprint is available 0|. 
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